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I. INTRODUCTION 



The static and dynamic chiral susceptibilities in the quark-gluon plasma are calculated within 
the lowest order perturbative QCD at finite temperature and the Hard Thermal Loop resummation 
technique using an effective quark propagator. After regularization of ultraviolet divergences, the 
£f*l • Hard Thermal Loop results are compared to QCD lattice simulations. 
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■ The main aim of rclativistic heavy- ion collisions is to unravel the basic properties of QCD associated with, e.g., its 
| phase diagram. Two of the main features characterizing the ground state of the theory, confinement and the sponta- 

■ neous breaking of chiral symmetry, are expected to cease in a possibly common phase transition at finite temperature 
^ ' and/or density [1]. The dynamical breaking of chiral symmetry is associated with the condensation of quark- antiquark 

ON . pairs in the QCD vacuum. As the temperature and/or baryon density increase, the QCD vacuum undergoes a phase 
transition to the chirally symmetric phase where the quark-antiquark condensate, the order parameter of the chiral 
phase transition, vanishes. The fluctuation of this order parameter is related to the associated susceptibilities [2]. 
The quantity of interest here is the chiral/scalar density susceptibility which measures the response of the chiral con- 
^vi | densate/scalar density to the variation of the current quark mass. The chiral susceptibility was recently measured in 

■ lattice QCD [1,3] with two light flavors. The chiral susceptibility has also been studied in chiral perturbation theory 
[4] , in the multiflavor Schwinger model for a small nonzero quark mass [4] , in the NJL model [5] , and using the Dyson- 
Schwinger equation [6]. In this paper, we will investigate the chiral/scalar density susceptibility following the same 
line as in the case of the quark number susceptibility [7] and the free energy of the quark-gluon plasma [8,9] within the 

Oh, Hard Thermal Loop (HTL) approximation [10], which selectively resums higher order corrections corresponding to 
2^ ■ medium effects such as screening, quasiparticles masses and Landau damping beyond the usual perturbation theory. 



II. DEFINITION 



> 

X 
b . 

A. General 

Let O a be a Heisenberg operator. In a static and uniform external field T a , the (induced) expectation value of the 
operator O a (0, x) is written [2] as 

Tr \0 CO x)e-^ n+n ^] 1 f 
cp a ee (O a (0, 50) = 1 ^H^+tt..)] J = y J ^ (O a (0, X )) , (1) 

where translational invariance is assumed and H ex is given by 

n ex = -J2 fd 3 xO a (0, x) T a . (2) 

The (static) susceptibility Xap is defined as 
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Xa0(T) 



d(j) a 



F=0 



= (3 J d s x(O a (0,x)Op(0,6)) , 



(3) 



assuming no broken symmetry (O a (0,x)) = (Op(0,0)} = 0. The quantity (O a (0,x)Op(0,())) is the two point 
correlation function with operators evaluated at equal times. 



B. Chiral Susceptibility 

The chiral susceptibility measures the response of the chiral condensate to the infinitesimal change of the current 
quark mass m + 5m. Here, T a corresponds to the current quark mass and the operator O a to qq. Then the static 
chiral susceptibility can be obtained [2] from 



Xc(T) 



d(qq) 



dm 



= f3 / d 3 x (q(0,x)q(0,x)q(0,d)q(0,d)) = (3 



I J d 3 x S(0, x) 



where S(Q, x) = ( ) is the static correlator of the scalar channel. The chiral condensate is given as 

Tr [qqe-^] _ 90 



(99) 



Tr \e-< m \ 



dm 



(4) 



(5) 



Here, H is the dynamical Hamiltonian of the system containing the quark mass, il = —(T/V) logZ is the thermody- 
namic potential, and Z is the partition function of a quark- antiquark gas. 

Taking the Fourier transform of the static correlator 5(0, x) = ( ), it can be shown that 



Xc(T) 



/ + OC 
-oo 



dw 
2^ 



S(w,p) , 



(6) 



where S(ui,p) is the Fourier transformed correlator. 

Applying the fluctuation-dissipation theorem [2] , one finds the static limit of the dynamic chiral susceptibility 



Xc (T) = lim (3 
p— »o 



du> 



■ Imll (u),p) 



, 2?r 1 - e~&" 

where H(uj 7 p) is the scalar self-energy. Using ImII(— co,p) = —lmU(uj,p), (7) can also be written as 

Xc(T) = -lim ^ / duj coth(^) Imn(w,p) . 

p^O 7T J Q \ 2 J 



(7) 



(8) 



III. PERTURBATION THEORY 



A. Static Susceptibility 

The chiral condensate (qq) at finite temperature can be obtained from the tadpole diagram in Fig.l as 



(qq) = -N c N f 



k 



d 3 k 



-Tr 



S(K) 



(9) 



where, Nf and A^ c are, respectively, the numbers of quark flavors and colors. S(K) is the free quark propagator given 
as 



S(K) = 



1 

— m 



(10) 
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Summing over the Matsubara frequencies fcoi we obtain 



(qq) = -2N C N f 



d 3 k 



(2tt) 3 E k 



1 - 2n F (E k 



(11) 



where n F {E k ) is the Fermi distribution function and Ek = Vk 2 + to 2 . 




Fig.l: Tadpole diagram determining the free static chiral susceptibility. 



Alternatively, the chiral condensate can be derived from the lowest order thermodynamic potential of the quark- 
antiquark system [11] 



n = —2NfN c T 



d 3 k 



(2 



log 



1 + e 



-0{E k -^ q ) 



+ log l + e - f){Ek+ ^ } 



(12) 



Using (5) the condensate (11) is reproduced. 

Now, the chiral susceptibility in the free case follows from (4) and (11) as 



X f c (T) 



d(qq) 



m=0 



2NfN c 



2N f N c 



d 3 k f P_ 



\-2n F (E k ) + 2/3-2 n F {E k ) 1 - n F (E k ) 



m=0 



d 3 k 1 

(2tt) 3 k 



1 - 2n F (k) 



(13) 



Note that, in contrast to the quark number susceptibility [7] where the zero temperature contribution vanishes due 
to the transversality properties of the vector channel [12], the static chiral susceptibility contains a quadratic ultra- 
violet divergence coming from the zero temperature contribution. Using dimensional regularization the temperature 
independent term in (13) disappears, as there is no scale associated with it, leading to 



xl(T) = - 



N f N c 
6 



(14) 



The next order contribution follows from the two-loop thermodynamic potential given in Ref. [11]. However, the 
second derivative of this expression with respect to m diverges at to = (see also Ref. [13]). The reason for this 
divergence can be seen in the following way: the second derivative with respect to the quark mass to of the two- 
loop thermodynamic potential, which is a quark loop with an internal gluon line, introduces two additional quark 
propagators. Hence the originating diagram looks like the polarization tensor containing an internal quark line. This 
diagram is known to be logarithmically infrared divergent in the case of a vanishing bare quark mass (see e.g. [14]). 
Hence, the static chiral susceptibility cannot be calculated consistently in usual perturbation theory beyond leading 
order but requires the HTL resummation. 



B. Dynamic Susceptibility in the Static Limit 



Next, we calculate the static limit of the dynamical chiral susceptibility from the lowest order self-energy diagram 
shown in Fig. 2: 
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U(P) = N f N c Tj2 



d 3 k 



S(K)S(Q) 



(15) 



Summing over fcoj we obtain 



Imn(w,p) = -NfN c ir 



d 3 k 



k ■ q + in 2 

E k E q 



k ■ q + m 2 

EkE q 



1 - n F (E k ) - n F {E q ) 5(co -E h - E q ) ~ 8{uo + E k + E q ) 



n F (E k ) - n F {E q ) 5{u - E k + E q ) - 5{u + E k - E q ) 



(16) 



K 



P 




P 



= P~K 



Fig. 2: Self-energy diagram determining the free dynamic chiral susceptibility. 



Combining (7) and (16) yields 



X f c (T) = 2N f N c 



2N f N c 



d 3 k 
(2tt) 3 



k 2 

1 - 2n F (E k ) + 2n F (E k ) 



TYi 

20 n F (E k ) 



1 - n F (E k ) 



m=0 



d 3 k r 



(2tt) 3 



1 - 2n F (k) + 2n%(k) 



(17) 



Obviously (13) and (17) do not agree. For, in general, the static susceptibility, given in (13) as the second derivative 
of the thermodynamical potential, does not coincide with the static limit of the dynamical susceptibility in (17) 
which satisfy the fluctuation-dissipation theorem. There is an exception: if the static susceptibility were associated 
with a conserved quantity such as the (quark) number density, the conservation law would ensure the equality of the 
two different quantities [2,15]. Thus, as shown in Refs. [2,7], the quark number susceptibility can be derived in two 
different ways as the imaginary part of the self-energy is proportional to S(u>) implying a conserved number density. 
However, the static chiral susceptibility is given by the derivative of the scalar/chiral density in (13) which obeys no 
conservation law, so it does not coincide with the static limit of the dynamic susceptibility in (17) as represented by 
the fluctuation-dissipation theorem since the imaginary part of the self energy in the static limit is not proportional 
to S(uj). 

Note that, in contrast to the static chiral susceptibility, there is no zero temperature contribution as there is an 
overall factor 0. The term containing no distribution function in (17) shows a cubic ultraviolet divergence. 

IV. HTL APPROACH 

Here, we want to calculate the chiral susceptibility beyond the free quark approximation by taking into account 
the in-medium properties of quarks in a QGP. In the weak coupling limit [g <C 1), a consistent method is to use the 
HTL-resummed quark propagator and HTL quark- meson vertex 1 if the quark momentum is soft (~ gT). However, 



1 In the case of a scalar meson there is no HTL contribution to the quark-meson vertex [16] 
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in the following we will use the HTL-resummed quark propagator also for hard momenta, in order to consider, at 
least to some extent, non-perturbativc features of the QGP such as effective quark masses and Landau damping. This 
strategy is in the same spirit as in the case of the quark number susceptibility [7] and the meson correlators [16], 
although this approach might not be consistent in the weak coupling limit [17]. However, since we do not want to 
restrict to the weak coupling limit, but consider temperatures close to T c , we do not aim at a complete leading order 
calculation, using here the HTL Green functions only for including medium effects in a gauge invariant way. 

Another, maybe more convincing, reason for studying the chiral susceptibility in the HTL approximation is that 
this approach has also been used for the free energy of the QGP, where different implementations of the HTL method 
and their validity have been discussed. The chiral susceptibility may serve as another test for these questions. 

In the rest frame of the medium the most general expression for the fermion self energy in the chiral limit of 
vanishing current mass has the form [18] 



Z(K) = a£-by , 



(18) 



where the scalar quantities a and b are functions of energy ko and the magnitude k of the three momentum. The 
HTL-rcsummed quark propagator, S*(K) = 1/[^T — T,(K)], can be expressed using the helicity representation in the 
massless case as [19] 



S* (k ,k) = 



k.j 



k.j 



2D+(K) 2D-(K) 



(19) 



with 



D±(k ,k) = (-k ± k) (1 + a) 



= -k ±k + 



in- 



TQi 



(20) 



where the thermal quark mass is given by m q — g(T)T/^/Q, and Q n (y) is the Legendre function of second kind. 
The zeros ui±(k) of D±(K) describe the two dispersion relations of collective quark modes in a thermal medium [19]. 
Furthermore, the HTL-rcsummed quark propagator acquires a cut contribution below the light cone (k 2 < k 2 ) as the 
HTL quark self-energy, contained in D±(K), has a non-vanishing imaginary part, which can be related to Landau 
damping for spacelike quark momenta resulting from interactions of the quarks with gluons of the thermal medium. 
For the coupling constant entering m q (T), a s = g 2 /(47r), we choose the two- loop result of the temperature dependent 
running coupling 2 [9] 



4-7T 



/3 2 L 



(21) 



where /3o = (11 — f-^/)) Pi = (51 — ^-Nf), L = log ^/z|/A|^.^ and ^ 4 <~ T is the renormalization scale for the running 
coupling constant. It should also be noted that the HTL-rcsummed propagator is chirally symmetric in spite of the 
appearance of an effective quark mass [18]. 

The HTL-spectral function of 1/D±(K), p± = — Im(l/D±)/ir, can be written as [19] 



p±{ko,k) = % 



2m 2 



5(k - u>±) + S(k Q + lo t ) 



+ P ± (k ,k)Q(k 2 -k 2 ) 



with 



(3±(k Q ,k) 



k(-k ±k) + m 2 q (±l-±^log*±£) 



1 2 




+ 


tt™2 ±k — k 
2 m q k 



2 ' 



(22) 



(23) 



where the first term describing the quasiparticle dispersion in (22) is due to the poles of the propagator and the second 
term containing Landau damping corresponds to the cut contribution for spacelike quark momenta. 



2 In this way, we have combined an infrared (HTL) and an ultraviolet (renormalization group) resummation scheme phe- 
nomenologically. A systematic combination of these techniques is not available yet. 
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A. Static Susceptibility 



Now, we compute the chiral condensate (qq) from the tadpole diagram in Fig.l, where we use the effective HTL 
quark propagator, yielding 



(qq) 



-N n Nf T 



fc 



d 3 fc 



.-Tr 



S*{K) 



(24) 



A similar calculation has been done in the case of the quark contribution to the free energy by Andersen et al. [9]. 

According to (4), but in contrast to the HTL approximation discussed above, we have to take into account a 
non-vanishing current quark mass. In the weak coupling limit (HTL approximation) there are two possibilities: 1) 
to > T, i.e., m 3> Shtl ~ gT and 2) to < gT. In the first case, Shtl can be neglected in the effective quark 
propagator compared to the hard scale to, implying S*(K) = S(K), which leads to the free susceptibility. In the 
second case, to can be neglected in the calculation of the HTL quark self energy, requiring hard loop momenta, i.e., 
£ ~ Shtl(^ = 0). Then, the effective HTL quark propagator becomes 



S*(K) = 



1 



${l + a) + b ln +m 



$(l + a) + &70 - m 



D 



m* 



with 



D(K) = K 2 {1 + af + 2fc 6(l + a) + b 2 = D + (K)D_(K), 



(25) 



(26) 



where D±(K) is given by the massless HTL-approximation in (20). Note that the current quark mass, m < gT, is of 
the same order as the HTL contributions a and b and cannot be neglected therefore. 
Summing over fcp, we obtain 



(qq) = 4mN f N c 



d 3 k 



dxp(x, k, to) [1 — uf(x)} 



where the spectral function corresponding to the massive HTL propagator is given by 

1 



p(x, k, to) 



"~ Im "n 

7r D — m z 



(27) 



(28) 



Following (4), the static chiral susceptibility in HTL-approximation is found as 
d(qq) 

d 3 k (uj 2 + (k)-k 2 



X h c (T) = 



dm 
ANfN c 



(2tt) 2 



2to2 



-ReDZ 1 (uJ+,k) [1 - 2n F (u)+)] 



(k) - k 2 



2m 2 q 



Rc J D+ 1 (w_,/c) [1 - 2n F {u 



+ 



dx [(3+(x, k)RcDZ 1 (x, k) + 0-(x, fc)Re£>- 1 (a;, k)] [1 - 2n F (ar)]| 



where 



ReD^^fco.fc) 



k 2 (-k ±k)+m 2 q k(±l- 



2k 



log 



fc+fcp 
k — ko 



) 



k{-ko ± k) + ml (±1 - ±^ log 







2 






+ 


k — ko 


\ 





Iml^^Qik 2 - k 2 ) 



(29) 



(30) 



In the infinite temperature limit, the first term in (29) reduces to the free static susceptibility given in (13) whereas 
the second and third terms vanish. The integrals without the distribution functions in the first and third term in (29) 
are ultraviolet divergent (quadratically and logarithmically) because of the asymptotic behavior of the quasiparticle 
dispersion relation u + {k), RcDjj 1 , and (3±. However, the second term is convergent since the plasmino dispersion 
relation, approaches k exponentially at large k. The quadratic divergence is temperature independent and can 

be removed by subtracting the zero temperature part, coinciding with dimensional rcgularization, or by subtracting the 
free chiral susceptibility (13) as discussed below. The remaining logarithmic divergence turns out to be temperature 
dependent as in the case of the gluonic part of the free energy [9]. This might be an artifact of the approximation 
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used here and might vanish if one goes beyond the one-loop HTL approximation. Here we will use dimensional 
regularization to remove this divergence following Ref. [9] . 

First we subtract the free susceptibility given in (13) from the HTL one in (29): 



S-rf = ^jf d**„, W + ^jf d,^^Re^(. + , fc)M . + ) 
_2^_/V £ r dfcfc 2 ^- - fc2 ) Re g-V-,fc) (1-2M«-)) 

TT 2 7 L 2m2 



+ 



47V N /*°° f 

^ dkk 2 dx [(3+{x,k)ReDl 1 (x,k)n F {x) + 0-(x,k)ReD+ 1 (x,k)n F (x)] 

n Jo Jo 

OAT AT f 00 

^ dkk 2 dx [p+{x,k)ReDZ 1 (x 1 k) + ^{x,k)ReDl 1 (x 1 k)] . (31) 

71 Jo Jo 

The integrals in 4 th and 6 th terms are, respectively, logarithmically ultraviolet divergent. Using a momentum cutoff 
k < A and energy cutoff x < A, the 4 th and 6 th terms can be seen to be proportional to m 2 log A. This is in contrast 
to the quark contribution to the free energy but similar to the case of the gluonic part of the free energy where these 
terms do not cancel [9] . Dimensional regularization will then replace the remaining logarithmic divergences by poles 
in d — 3, which finally will require a counter term for cancellation [9]. 

As pointed out in Ref. [20] it is important to use also ci-dimensional expressions for the HTL propagator. Here 
we will follow Ref. [21,22] and expand the chiral susceptibility up to order trig. The angular average leading to the 
Legendre functions in (20) can be generalised to d = 3 — 2e spatial dimensions in Minkowski space as 



where the weight function w(e) is 



r(§)r(i- e ) 



Eq. (32) can also be written as [21] 



Jo 



k 2 I k 2 



T d (K)=w(e)l dc (l_ c ^ r ^_iL_ = ^_ : iL_^. (34 ) 
Now, it is convenient to define the functions Aq(K) and As(K) as 



A = k - ^%(K) , 



Tfl 

A s = k + -jl{l-T d (K)), (35) 

which are related to the HTL-resummed quark propagator (20) by 

D±(ko,k) = -(A TA s ). (36) 
Using (4), (24), and (25) the chiral susceptibility in one- loop HTL order reads 

X h c (T) = -AN f N c ^ — L_, (37) 



where [22] 



{*} 
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with the arbitrary renormalization scale /i of the (MS) renormalization scheme and Eulcr's constant 7. 
Expanding Xc(T) hi a power series of m 2 /T 2 , and keeping terms up to second order, we obtain 



X h c (T) = -AN f N; 



1 

k 2 



2ml 



1 



K 6 k 2 K 4 



1 %(K)+ : 



k 2 K 4 



k 2 k 2 K 2 



(39) 



Using the sum-integrals evaluated in Appendix A, we end up with 



X h c (T) 



NfN Crr2 



1 + 



-( 

IT 2 \ 



AttT 



- + 2 7 + 41og2 -f + 



21C(3) 



T 2 4tt 4 



7r 2 \ mf 



5 + - log2 - — ) ^ 



(40) 



The chiral susceptibility is logarithmic divergent at order m 2 , but finite at order m q . We choose a counter term within 
the minimal subtraction renormalization scheme just to cancel the pole in 1/e: 



N f N, 



f ly c 2 

m„. 



2tt 2 



(41) 



The complete expression for chiral susceptibility is then given by 



X h c (T)=x f c (T) 



1 + ^( log 4£r + ^ + 21og2 ) 



Q 
y2 



21C(3) 

4 7T 4 



5 + ylog2 



7T \ m , 

rp4 



(42) 
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Fig 3: The HTL static susceptibility (free susceptibility subtracted) as a function of temperature for Aj^g = 300 MeV, 
and Nf = 2 with the choices of the renormalization scale fi = 2nT and AitT for up to order m 2 (left) and up to order 
mi (right). 



Note that the m 2 -term is the first finite a s -correction to the free chiral susceptibility, since the two-loop contribution 
in usual perturbative QCD diverges for m = 0. 

The quantity \x^{T) — x£(T)|/T 2 as a function of T is shown in Fig. 3 for two different choices of /z, where we 
have also chosen /i4 = in a s entering m q . Obviously the nig contribution is only a small correction. Whereas in 
the m 2 contribution only the effective quark mass (pole contribution) enters, there is a Landau damping contribution 
coming from 7^ in the m q term (see (39)). To understand the behavior of the susceptibility in Fig. 3 we ascertain 



that the quantity 



IXc 



xl\/T , plotted in Fig. 3, depends on the temperature only via the temperature dependent 



running coupling constant, i.e. \Xc — Xh\/T — F\g(T)]- This can be seen by scaling the variables k and x in the 
integrals of (29) by m q . It holds also after dimensional rcgularization, if a renormalization scale proportional to the 
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temperature is adopted in the latter case (see e.g. (42)). This dependence has also been verified numerically, where 
a constant, temperature independent chiral susceptibility has been observed, if a constant, temperature independent 
coupling has been employed. In other words, the decrease of the susceptibility with temperature in Fig. 3 is entirely 
caused by the temperature dependence of the coupling constant. Due to asymptotic freedom, g(T — ► oo) — > 0, the 
HTL susceptibility reduces to the free one in the infinite temperature limit, i.e. F[g — 0] = 0. 

In lattice simulations a peak around the critical temperature associated with chiral symmetry restoration has 
been observed. In contrast to lattice QCD the HTL approximation, although taking into account medium effects 
in the plasma due to interactions, does not contain any physics related to the chiral restoration, which is a truly 
non-perturbative effect. In Fig. 3 we observe a strong increase of the magnitude of the HTL chiral susceptibility 
towards low temperatures, similar as in lattice simulations above T c . However, this increase is caused only by the 
fact that for our choice of the renormalization constant n the HTL result is above the free one, \x^\ > |x£|, and 
that the temperature dependent coupling increases at smaller temperatures. If we choose a small enough value for 
Mi IXcl < IXcl holds according to (42), and the HTL susceptibility |xjf| approaches the free susceptibility |x£| from 
below at large temperatures. This indicates the strong sensitivity to the choice of the regularization procedure and the 
renormalization constant similar as in the case of the free energy [21,22]. Hence we understand the difference between 
the lattice and the HTL approach qualitatively but the different ultraviolet behavior of the HTL and the lattice 
result, which does not contain a temperature dependent logarithmic divergence [23], renders a direct, quantitative 
comparison of the both approaches difficult. 



B. Dynamic Susceptibility in Static Limit 



Here, we need to calculate the imaginary part of the self energy diagram given in Fig. 2, where both bare quark 
propagators are replaced by effective HTL quark propagators. It is sufficient to restrict to the massless case, i.e., (19). 

Apart from constant factors the dynamic chiral susceptibility agrees with the scalar meson correlator at r = (i 
[compare (2.3) of Ref. [16] with (8)]. Hence, we will not repeat the derivation here but just present the final results. 

As in the case of the meson correlator [16], the dynamic susceptibility can be decomposed into pole-pole, pole-cut 
and cut-cut contributions, 



Xc(T) = X P c P (T) + X P c c (T) + X C c(T). 



(43) 



The pole-pole contribution reads 

f d 3 / 



X PP (T) 



3 k_ 

? 



4m* 



1 — 2n F (uu + ) + 2n F (w + ) 



+ - 



.(fc)-fc 2 



1 - 2n F (w_) + 2n 2 F (u-) 



+ 



[<4(fc) - k 2 }[u 2 _(k) - k 2 } 



2m* 



n F (u) + )[l — njr(w_)] + nj?(w_)[l — n F (io + ] 



']}■ 



(44) 



In the infinite temperature limit the plasmino mode decouples from the medium, causing the second and third terms 
to vanish, whereas the first term reduces to the free susceptibility given in (17). The pole-cut contribution is given as 



XT{T) 



2N f N c p 



\ Jo 



duu coth 



//3w\ f d 3 fc 
I 2 J J (2tt)3 



m' q Jo \ 2 J J V*) 

{0[fc 2 — (lu — LU + ) 2 }riF(uJ — uJ+)riF(u + )f3 + (uj — uj + , fc)(w"j_ — k 2 ) 
+@[k 2 - (lu- LU-) 2 ]n F {uj - w_)n F (w_)/3_(o) - w_,fc)(w 2 - k 2 ) 
+<d[k 2 - (lu + LU-) 2 ]n F (LU +u>-)[l - ti f (lu-)](3 + (lu + w_, k)(LU 2 _ - k 2 ) 
+6[k 2 -(lu + LU + ) 2 ]n F (LU + lu + )[1 - n F (uj + )]f3-(LU + lu + , k)(uj 2 + - k 2 )} , 

and the cut-cut contribution as 

Xf{T) = 2N f N c P J °° cLu coth (^j J ^3 j^dx n F (x)n F (Lu - x)Q(k 2 - (x - lu) 2 ) 

P+(x, k)fi + (uu — x, k) + (3-(x, k)/3-(u) — x, k) . 



(45) 



(46) 
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The dynamic chiral susceptibility in the 1-Loop HTL approximation, considered here, has a cubic ultraviolet 
divergence, which can be removed, for example, by subtracting the free result (17). We will not discuss the dynamic 
susceptibility further on as it is not directly related to QCD lattice results [3]. We will just mention here that the 
(ultraviolet divergent) dynamic susceptibility, (44) to (46), normalized to the free case (17) is equal to unity, as can 
be seen from the scalar meson correlator at r = f3 discussed in Ref. [16]. 

V. CONCLUSIONS 

In the present paper we have considered the chiral susceptibility of the quark-gluon plasma within lowest order 
perturbation theory (free case) and the HTL approximation. We have shown that there are two possible definitions 
of the chiral susceptibility leading to different results due to a missing conservation law associated with the chiral 
susceptibility in contrast to the quark number susceptibility [7]. Using the definition based on the derivative of the 
chiral condensate with respect to the current quark mass, we have computed this static chiral susceptibility within 
perturbative QCD and using an effective HTL rcsummed quark propagator. In the case of the usual perturbation 
theory we found that the chiral susceptibility cannot be computed beyond the leading order (free case) as it diverges 
due to setting the bare quark mass equal to zero at the end. Similar as in the case of the free energy [9] but in contrast 
to the quark number susceptibility [7], the chiral susceptibility using the HTL approach is ultraviolet divergent. Even 
after subtracting the free susceptibility or the zero temperature part, a logarithmic temperature dependent singularity 
remains as in the case of the gluonic part of the free energy [9] . Applying dimensional regularization the HTL chiral 
susceptibility can be larger or smaller than the free one depending on the choice of the renormalization constant /z. 
Here we expanded the chiral susceptibility up to order m A q . The m 2 contribution is the lowest order correction to the 
free susceptibility, which has nom = divergence. The temperature dependence of the HTL chiral susceptibility is 
entirely determined by the temperature dependence of the running coupling constant that enters through the medium 
effects (quasiparticle quark mass, Landau damping) contained in the HTL approximation. No indication of a peak at 
the critical temperature is found, as seen in lattice QCD simulations [3] . This indicates that this peak is due to truly 
non-perturbative effects such as chiral restoration, not included in the HTL approach. 
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APPENDIX A 



1. One-loop sum-integrals 



The one-loop fermionic sum-integrals necessary for our purpose in (39) are: 

,C'(2) 




^ k 4 K 2 ^ k e K 2 



1+ (2 7 + 21og7r-2^)e + 0(e 2 ) 



L+(2 + |log2 + 2^), hCV) 



(Al) 



(A.2) 



tt) 2 \4ttTJ 



(K 2 ) 2 (4tt) 



2c 



l - + 2( 7 + 2 log 2) + f + 4 7 log 2 + 4 log 2 2 + 4 7l ) , + 0\ * ; ) 



(A.3) 



k 2 K A UttTV 2(2tt) 4 j 



2(27r) 4 T 2 



I f (4+^log2 + 2ffl) ( fCV) 



(A.4) 
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21C(3) 



(K 2 ) 3 128tt 4 T 2 \4ttT 



(A.5) 



All these above sum-integrals are calculated using standard contour integration prescription. The number 71 in 
(A. 3) is the first Stieltjes gamma constant. 



2. One-loop HTL sum-integrals 



In this subsection we will compute the one-loop sum-integrals in (39) involving the HTL function Td(K) using the 
method developed in Ref. [22]. 

The following sum-integral in (39) can be decomposed using (34) as 



k 2 K 4 



UK) 



1 -c 2 



c 2 - c 3+2e 



k 2 K 4 \ (1-c 2 ) 2 I c ^ } k 4 K 2 ' 



(A.6) 



The one-loop sum-integrals are given in (A. 2) and (A. 4). The angular averages are performed using (34), leading to 

1 \ 1 



1-c 2 



2e 



1, 



(A.7) 



C — C T 



1 

~4~e 



log 2 



3 it 2 



2 log 2 -log 2 2 e + 0(e 2 ). 



(1-c 2 ) 2 / c 4e ' \4 "° J ' V 4 6 

Once these expressions are combined, the pole in 1/e cancels and we obtain a finite result: 



1 7<(3 » 



k 2 K 4 



4(2tt) 4 T 2 



Following the same strategy, the remaining sum- integral in (39) involving T 2 (K) can be written as 



i k 2 k 2 K 2 



-cf + c 2 (-l + c|) + cf - c 2 (-l + c|) 
(c 2 -c 2 )(-l + c 2 )(-l + c 2 ) 



k 2 
ft o 

fc 6 if 2 



where e has been set to zero, as this integral is finite. The result is given as 



1 tfm- rm 



^ } k 2 k 2 K 2 dK ' (2^) 4 T 2 



21og2 

B 12 



(A.8) 



(A.9) 



(A.10) 



(A.11) 
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